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ABSTRACT 
A method is developed to analyze and compare the effectiveness of 
ASW vehicles. The measure of effectiveness is the probability that the 
vehicle, after detecting a submarine with passive sensors, can transit 
to the contact area and re-establish contact with the submarine. A com- 
puter simulation is developed and an example using three hypothetical 


ASW vehicles is illustrated. 
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|i Introduction. 

As a result of increased emphasis currently being placed on ASW 
weapon systems, it is advantageous to possess a computer simulation 
ei of analyzing and comparing proposed ASW vehicles and opera- 
tional ASW vehicles. The computer simulation for the comparison of ASW 
vehicles developed in this paper is intended to be used for this pur- 
pose. 

The problem at hand is to develop a method for comparing ASW vehi- 
cles under operational conditions using vehicle, submarine, and tactical 
situation parameters. The simulation allows one to compute a proba- 
bility which is used as a measure of effectiveness. This probability, 
hereafter referred to as the "probability of success,"' is defined in 
this paper as the probability that an ASW vehicle successfully relocates 
a submarine at an estimated position given that a detection has occurred. 

A description of the simulation and its limitations is made in Sec- 
tion 2. It is very important to understand the simulation and the 
background assumptions before the results of the simulation are ana- 
lyzed. After the problem is developed mathematically in Section 2, a 
discussion and an illustration is given in Section 3 to demonstrate how 
to use the computer program. In Section 4, several applications are 
discussed. Finally a comparison is made of three hypothetical ASW 
vehicles--a High Speed Surface Vehicle, a VERTOL Vehicle, and a Seaplane- 
Type Vehicle--using the results generated from the computer simulation. 
The output of one computer run is illustrated in Figures 5 and 6. Curves 
for the probability of success versus range to the submarine for each 
vehicle are displayed in Figures 7, 8, and 9. These curves are used to 


analyze and compare the systems. 








The results of these comparisons illustrate the usefulness of the 
simulation in comparing the performance of different vehicles and for 
studying the effect of various parameters such as expected sensor range 


on the probability of success. 


Ze Simulation Description and Mathematical Development. 

Effort has been directed toward the development of a computer simu- 
lation that represents a real world conflict between an ASW vehicle and 
a submarine. The setting is one of an ASW vehicle assigned to a patrol, 
surveillance, or barrier mission. Initial detection is assumed to have 
been made using passive sensors. The vehicle then proceeds as rapidly 
as possible to an estimated position (EP) in order to relocate the 
submarine. The EP is the predicted position of the submarine when the 
actual position is not known, but the vehicle's sensors indicate the 
presence of a submarine. The objective of the vehicle is to compute an 
EP of the submarine, close the position at maximum speed, and attempt 
to relocate the submarine using active sensors. No information about 
the actions of the submarine is available while the vehicle closes to 
the EP. 

Certain other restrictions are assumed in the model. Ahead thrown 
weapons and nuclear weapons are not available for use by the vehicle. 
Multiple targets are excluded. A passive detection range of at least 
five nautical miles is necessary; otherwise, an active sensor would be 
used and the problem of transiting to an estimated position is irrele- 
vant. Also, the vehicle is assumed to lie motionless in the water until 
detection has been achieved to insure passive sensor capability. 

The simulation is developed for two different tactical situations. 
The first case is that of a submarine which is unaware of the activity 
of the vehicle. The second case is that of a submarine which is alerted 
by the activity of the vehicle. Common to both of these cases is a 


maneuvering board problem which will be discussed first. 








The maneuvering board solution consists of describing geometrically 
the positions and motions of a vehicle and a submarine with time. The 
passive sensor of the vehicle provides the last known position, course 
and speed of the submarine, and the bearing and range to the submarine 
from the vehicle. Knowing these values and the vehicle speed, an EP, 
transit time of the vehicle, relative closing speed, and the distance 
traveled to EP are computed. For overall ease in reading the same nota- 
tion is used in the mathematical development that has been used in the 
computer simulation (Appendix C). 

The solution is computed by solving a velocity vector problem. The 
vehicle is located at the origin (0,0) of a rectangular coordinate sys- 
tem. One side of the velocity triangle (Figure 1) represents the 
estimated velocity vector of the submarine (U1) which is re-established 
at the origin. The velocity vector for the vehicle (Z) is unknown in 
direction, but known in magnitude. The remaining relative velocity vec- 
tor (RELSPD) represents the relative course and speed of the vehicle 
closing the submarine. Only the direction of RELSPD is known which is 
parallel to the true bearing line from the vehicle at (0,0) to the sub- 
marine at (X2, Y2). Then the intersection (X3, Y3) between the relative 
velocity vector and the vehicle velocity circle is computed. 

In order to solve the maneuvering board problem various points and 
vectors in the rectangular coordinate system must be computed. The 
values for the components of Ul are: 


UL x SIN (THETA) (1) 


Xl 


YL = UL x COS (THETA) . (2) 


The Last known position (X2,Y2) of the submarine is given by: 
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Figure \. “The Maneuvering Roara Ssiotiar 








X2 = RAI x SIN (BETA) (3) 


Y2 


RAI x COS (BETA) (4) 
where BETA is the bearing in degrees of the submarine from the vehicle 
and RAI is the range in nautical miles of the submarine from the vehicle. 
The RELSPD vector is represented by the line through the point (Xl, Yl) 
and parallel to the line connecting the origin and (X2, Y2). The equa- 


tion of this line is: 


Y3 = = (X3-X1) + YL. (5) 


The RELSPD vector intersects the vehicle velocity circle at the point 
(X3, Y3). The equation of the vehicle velocity circle is: 

¥32 + x3? = y2_ ) 
where V is the average vehicle transit speed. The point (X3, Y3) is de- 
termined by solving equations (5) and (6) simultaneously. The direction 
of the vector from the origin to (X3, Y3) represents the course of the 
vehicle when traveling inbound to EP. The relative speed vector is rep- 


resented by the vector between (Xl, Yl) and (X3, Y3) and its magnitude is; 


RELSPD = 4] (X3-X1)* + (Y3-Y1)2. (7) 


The vehicle's transit time (TD) to EP is defined as: 
TD = RAI/RELSPD , (8) 
Let the vehicle reaction time (RT) be the following sum: 


RT = PTOT + TOT + CLOT + RELT ; (9) 
where , PTOT - pretake-off time which commences when the vehicle can 
no longer observe the submarine and includes decision 

time, warm-up time, and sensor retrieval time, 


TOT - take-off time which includes take-off roll and time to 
take-off into the wind, 





CLOT - land time which includes time to land and taxi time 
needed to arrive at EP in order to deploy active sensors, 


RELT - target relocation time which includes the time necessary 
to deploy sensors and ends when the submarine is re- 
located. 

Any of these times can have a zero value. For example, a helicopter 
could have TOT equal to zero if it needs no take-off time. Also, the 
times are mutually exclusive so that no time interval is added twice. 
Obviously the blind time (TB) is: 

TB = RT + TD. (10) 

Evasion for the alerted submarine begins when blind time commences. At 
that time the submarine is assumed to be at (X2, Y2). If RT is equal to 
zero, the submarine is at (X2, Y2) when evasion begins; but if RT is 

not equal to zero, the submarine is not actually at (X2, Y2) when eva- 
sion begins because it takes RT for the vehicle to commence moving 
towards EP. During the reaction time the vehicle is motionless, but 

the submarine is moving away from (X2, Y2). Because of the small dis- 
placement of the submarine from (X2, Y2) due to RT, the assumption that 
the submarine is at (X2, Y2) when evasion begins does not significantly 
affect the results of the simulation. 

It is necessary to determine the estimated position (EP) of the 
submarine after a time interval of TB and a submarine speed U3. The 
equation of the estimated distance (DIST) traveled by the submarine is; 

DIST = U3 x TB. (11) 


The EP is also the position at which the vehicle estimates that an in- 
tercept will be made with the submarine. The components of EP relative 
to the origin are: 


X4 = DIST x SIN (THETA) + X2 (12) 





Y4 DIST x (COS (THETA) 2 (13) 


The equation for the distance DISDAT) traveled by the vehicle to EP is 


DISDAT = 1 xo ya? (14) 


The values for TB and DIST are very important in determining the proba- 


bilities of success for the two cases developed below. 


Case I. The unalerted submarine; 

It is conceivable that the vehicle could transit to EP without the 
submarine being cognizant of the impending danger. In this case the un- 
alerted submarine does not evade and proceeds at cruise course and veloc- 
ity. Prediction errors for the submarine course and velocity are made 
which are characteristic of the vehicle's equipment. The errors are 
assumed to be distributed normally about the predicted course and veloc- 
ity. Because of long blind times expected, the predicted positions of 
the submarine are not accurately represented by a circular-normal distri- 
bution, and a numerical integration must be performed. Probabilities of 
success are computed for each of a representative distribution of pre- 
dicted positions consisting of 25 points or cells (see Figure 2). This 
concept has been adopted from fails 

The last known submarine position (X2, Y2) is established at the 
Origin of a rectangular coordinate system. Since DIST is the distance 
traveled by the submarine during TB, then EP is (DIST,O) in this coor- 
dinate system. 

Using the prediction errors, a submarine position is computed for 
each of the 25 cell midpoints. The distance (VV) and angle (THETER) for 


each cell are computed using the submarine cruise speed (U1), the submarine 
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Figure 2. Numerical Integration For The Unalerted Case 








speed error (UERR), the submarine course error (THDEL), and the blind 
time (TB). The equations for VV and THETER are: 
VV = (U1-AA(IL) x UERR) TB (15) 
THETER = AA (J) x THDEL , (16) 
with AA(I) and AA(J) = 0, £ .52, + 1.28 depending on which cell is being 
considered. These numbers can be obtained from any normal density tables 
[2]. 


The equations for the submarine position (SUBX (I,J), SUBY(I,J)) at 


3 


the center of each of the 25 cells are; 
SUBX (I,J) = VV x COS (THETER) (17) 
SUBY (I,J) = VV x SIN (THETER) . (18) 
A bias is computed between EP and the center of each cell. The equation 


for the bias of the (i,j) cell is: 


SBIAS (I,J) = | (suex(1,3)-prsT)? + (SUBY(1,J))2 . (19) 


Various errors have to be considered to give the simulation realism. 
It is assumed that the navigation error and sensor locating error are 
additive about the point EP. The distribution of the vehicle location 
about EP due to these errors is assumed to be circular normal. 

The standard deviation of navigation error (SIGNAV) is: 

SIGNAV = PRCNIN x RAI | (20) 
where PRCNTN is a percentage of the last known range (RAI). The stand- 
ard deviation of sensor locating error consists of the standard devia- 
tion of bearing ( SIGSBR ) and the standard deviation of range (SIGSR) 
which are defined by: 

SIGSR = PRCNTS x RAI (21) 


SIGSBR = RAI x SIN (BRGER) (22) 


10 








where BRGER is the bearing error and PRCNTS is the percentage error of 
the last known range. Since (21) is greater than (22), because of the 
passive sensor, an elliptical normal distribution exists; but it is 
approximated by a circular normal distribution for ease of computation, 
This is accomplished by equating the area of an ellipse to the area of a 


circle. The result is an approximate standard deviation (SIGAPR). 


= 2 » with 
TT Ola a 2! ae 


YE = SIGSBR (23) 


OYE = SIGSR 


.°. sicapr = Ocrr =7/Oxe OVE (24) 


Therefore, the total standard deviation of error (SIGJT) is the sum 
of the squares of the approximate standard deviation of error and the 


standard deviation of navigation error according to the equation; 


SIGJT “7 (SIGNAV)2-+ (SIGAPR)2 . (25) 


Other error models were considered in order to more accurately represent 
the errors involved. Due to lack of time the error model described 
above, although not the best, was adopted. 

To compute the probability of success for the unalerted submarine 
the circular coverage function is used, [1]. In order to use this funce- 
tion, the parameters SBIAS(I,J), SIGJT, and EXPRNG are required. The 
first two have been defined above. The term EXPRNG, which is defined as 
the expected range of the active sensor used by the vehicle, is a charac- 


teristic of the sensor employed. The quantities SBIAS(I,J) and 


SIGJT 


EXPRNG 


SIGIT are used in the circular coverage function to compute the proba- 


bility of success (PK(I,J)) that a circular disk of radius EXPRNG will 


ll 








cover a point (SUBX(1,J), SUBY(I,J) from the EP since the probable posi- 
tion of the disk is described by a Gaussian distribution. This proba- 
bility PK(I,J) is the probability of success for the ijth cell. Since 
it is equally likely that the submarine is in each of the 25 cells, :.': 
the probability that it is within each cell is .04. The probabilities 
of success for each cell are summed and multiplied by .04 to give the 
"probability of success (PROB) against an unalerted submarine." The 


probability of success is: 


PROB = «| (PK(I,J) x .04) (26) 
J 


I 


Case II. The alerted submarine: 

The submarine is aware of the presence of an ASW vehicle and evades 
by changing course and speed. The speed changes are characteristics of 
the submarine but the course changes are determined in number and degree 
by the situation being simulated. A probability of success is computed 
for each evasion turn and an arithmetic average taken over the number 
of evasion tactics used. 

This case is designed to give practically full evasion to the sub- 
marine (see Figure 3). The submarine evasion turns are limited to a 
maximum turn of 90 degrees to either side of the estimated submarine 
course. No additional turns are granted after the initial turn is exe- 
cuted. For turns greater than 90 degrees the submarine would deviate 
from its initial track to such an extent that it would not be able to 
make up the lost time necessary to make good a mission speed of advance. 

The range of values for the set of evasion turns is determined by 


the situation being simulated. For example, a slow submarine would 
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almost always turn away from the vehicle when the vehicle's presence is 
known. There are two limiting turns which bound the set of evasion turns 
to the right and the left. The bounds must be no greater than 90 degrees 
to either side of the estimated submarine course. PSI and DELTA are the 
left bound and right bound, respectively. The range of values for the 
set of evasion turns is determined by (PSI-DELTA). The interval (ETA) 
between evasion turns is ETA a » NU being the number of eva- 


sion turns considered. If the set of turns encompasses turns to the 


right of the submarine course then ETA ( - In Figure 3, 


PSI = 060 degrees, DELTA = 330, and ETA = 060 (360-330) = 15°, 
6 


Therefore, the evasion turns considered for this example are 330°, 345°, 
360°, 015°, 030°, 045°, and 060°. 

In the alerted submarine case the submarine course and velocity pre- 
diction errors that were present in the unalerted case are ignored. 

Since in this case a set of evasion turns are selected which are only sub- 
jective predictions, it is valid to exclude prediction errors because 

they are comparatively insignificant to the computations. However, ve- 
hicle sensor and navigation errors that were assumed in the unalerted 
case still prevail. 

Initially, the submarine is at cruise speed. When evasion commences, 
the assumed submarine speed is increased to fast speed. The character- 
istics of the submarine determine how long the submarine remains at fast 
speed (length of time is Tl) after which it slows to silent speed. The 
time intervals of acceleration and deceleration are considered negligible 
when compared with the long blind times encountered. 


The alerted case uses a rectangular coordinate system with the last 
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known position of the submarine (X2,Y2) at the point of origin when eva- 
sion begins. From this point positions (X(N), Y(N)), which are the 
positions at the end of each of the N evasion turns, are calculated. 
(See Figure 3.) The time to turn (TIT), evasion turn in degrees (AB), 
and the turning radius (TR) are needed to compute X(N) and Y(N). The 


equations for TIT, X(N), and Y(N) are computed as follows: 


TTT = are length _ TRx AB 
sub speed U3 


If TI ZS TB, then 
X(N) = TR x SIN(AB)+U3 x (TB-TTT) COS(AB) (27) 
Y(N) = TR- TR xCOS(AB)+U3 x (TB-TTT) SIN(AB) . (28) 
If TI < TB, then 
X(N) = TRx SIN(ABH-U3x (TB-TTT) COS(AB)+ (TB-T1)U2xCOS(AB) (29) 


Y(N) = TR-TRxCOS(AB)+ U3 x (TB=-TTT)SIN(AB) + (TB-T1L) 
U2x SIN(AB), (30) 


The time Tl is the length of time the submarine remains at fast speed (U3); 
after which it changes to silent speed (U2). Equations (29) and (30) 
account for the speed change occurring before TB has terminated. The 
point (X(N),Y(N)) can be considered the submarine's actual position after 


th evasion turn is taken. 


the N 
The estimated position EP is the predicted position of the submarine 

after the interval of blind time has elapsed. Therefore, EP is (DIST,O). 

The distance between each (X(N),Y(N)) and EP constitutes a bias (BIAS(N)). 


This is an important factor in computing the probability of success. The 


bias for the NCP evasion turn is: 


BIAS(N) = ay (X(N)=-DIST)* + (Y(N) )2 : (31) 


As in the unalerted case, the parameters BIAS, EXPRNG, and SIGJT must be 


15 








computed (EXPRNG and SIGJT are the same as in the unalerted case) in 
order to use the circular coverage function. The circular coverage func- 
tion is used to compute the probability (PPK(N)), for each N, that the 
vehicle can regain contact at EP given the total standard deviation of 
error (SIGJT), the bias (BIAS(N)) between the terminal point of the nth 
evasion and EP, and the expected active sonar range (EXPRNG). 

The two ratios of the parameters that are needed for the circular 


coverage function are BIAS(N) and EXPRNG | with these ratios the proba- 
: SIGJT SIGJT 


bilities, PPK(N), are computed for each N, and then all PPK(N), N= l, 
~ee, NU+1, are summed and averaged over all NU+1. The average proba- 
bility is: 

NU+ 1 

PAVG = y = Dopey (32) 
— NUtlL 
It might be appropriate to use the option to weight each of the eva- 

sion turns. If it is believed that the turns within the set of evasion 
turns are not Seat likely, then the probability vector UNU(N), N#l, 
-2o, NU+1 can be used which gives probabilities that each of the N 
evasions is actually taken. If they are not equally likely, then the 
“weighted probability of success (WTPAVG) against an alerted submarine" 


is computed. The equation for WTPAVG is 


NU+1 


(PPK(N) x UNU(N)) 
WIPAVG = ¥- (33) 


(NU+1) 


The results of the mathematical development have produced methods 


of computing the probability of success (PAVG) for an alerted submarine 
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with equally Likely evasion turns, the probability of success (WTPAVG) 
for an alerted submarine with evasion turns that are not equally likely, 
and finally the probability of success (PROB) for an unalerted submarine. 


These probabilities are used to analyze and compare ASW vehicles. 
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a Input Data Rules, 

A complete run for the computer program (Appendix C) consists of 
25 input cards. Each input data card represents the value of one vari- 
able except for the last two cards, one of which contains the values for 
the AA(i), i#l, ...,5 and the other which contains the UNU(N) values, 
N=l, ..., NU+1. 

One of the first functions of the program is to set the indices 
NO=1 and M=l. (See Appendix B--Flow Charts.) The index NO is used to 
number the pages of the program output. The index M is compared with 
NEXT (initially set equal to one). The comparison dictates whether the 
program will continue or terminate. The value for NEXT determines the 
number of runs to be executed and is an input to the program. The compu- 
ter then reads the input data cards from SUBROUTINE INPUT after which 
the 25 values are printed on page 1 of the output. The index NO is in- 
creased by one. The programmed computations are performed and the 
results are printed on page 2 of the output (Figure 6). The indices NO 
and M are then increased by one. The index “ee compared to NEXT and 
if NEXT is greater than M, SUBROUTINE INPUT searches for more input data 
cards. Since the program is used to compare and evaluate systems, only 
a few variables will be changed for each successive run. However, more 
or all the variables can be changed as long as the input procedure is 
followed. For example, to compare the three vehicles in Table 1, first 
place the parameters for the High Speed Surface Vehicle, RAI = 5.0, and 
a value for NEXT in the data input for Run 1. (See Figure 4.) The next 
run consists of a Locator Card which initiates the reading of Variable 


Change Card, which in turn changes the value of RAI to 10. The subsequent 
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card would be the Information Card, which indicates no further changes 
for the second run -- start computations. For this run the variable RAI 
has changed and all other parameter values remain the same. 

If another variable change were desired, another Locator Card, 
Variable Change Card, and Information Card would be placed after the 
RAI Variable Change Card as indicated in Run 3 in Figure 4. 

A run following the initial run takes a minimum of three cards. If 
only one variable is being changed (e.g. RAI = 10), then three cards are 
necessary. This procedure is followed for RAI = 5, 10, ..., 60 for each 
system in Table 1 in order to compute probabilities of success for each 


vehicle over the entire RAI range. 
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4. An Example and Some Applications. 

The simulation possesses the capability of being used as a tool to 
analyze and compare ASW vehicles as explained in the previous sections. 
In order to illustrate this capability, the parameters Listed for each 
system in Table 1, and ranges RAI=5, 10, 15, ..., 60 nautical miles 
were used as input data in the computer program (Appendix B) and pro- 
cessed in the CDC-1604 computer. The program is written in FORTRAN-60. 
The flow chart is illustrated in Appendix C. 

A sample output for the Seaplane Type Vehicle using the parameters 
of Table 1 and RAI#45 is shown in Figures 5 and 6. The former is a list 
of all the parameter values for the indicated run. This list enables 
one to insure that the correct parameters for that run are being properly 
read into the Senter: The latter page contains the probabilities of 
success for one run. The probabilities for each RAI are plotted on 
graphs. Each vehicle is depicted on an individual graph. Each graph 
has six curves -- three drawn from values calculated from using a slow 
submarine as an adversary and three drawn from using a fast submarine as 
an adversary (see Figures 7, 8, 9). 

The parameter values in this sample are strictly hypothetical and 
demonstrate only the use of the simulation. However, with these chosen 
parameter values it can be seen in Figures 7, 8, and 9 that throughout 
the range of RAI the High Speed Surface Vehicle has probabilities of 
success considerably lower than the other systems. The slower the speed 
of the vehicle, the lower the probability because of a low closing rate 
and long blind time. The probabilities of success are higher at shorter 
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but a crossover point occurs at approximately RAI=25, where the Seaplane- 
Type Vehicle has higher probabilities of success. The crossover is a 
result of the values selected in the sample for the vehicle speed, reac- 
tion time, and expected sonar range. From these results the best vehicle 
can be chosen and the other vehicles investigated for improvements. 

In addition to the example stated above several applications are 
suggested for the simulation. One application is to use the simulation 
to determine what parameters or characteristics the vehicle and the 
sensors must have to obtain a given probability of success. For example, 
if the state-of-the-art for active sensors dictates that five nautical 
miles is the maximum range expected for active sensors for the next ten 
years, then it would be advantageous to find values of vehicle speed, 
reaction time, navigation errors, and prediction errors that can be tol- 
erated for the vehicle and still maintain a given probability of success. 
The development of a new ASW vehicle could be directed towards achieving 
these parameter values. 

Another application can be made if a new sensor is being proposed 
for a given vehicle. The values for the expected active sensor range 
and the sensor errors can be computed for a given probability of success. 
These parameter values can be used as design limits in the development 
of the proposed sensor. 

Finally the simulation can be used if intelligence indicated that 
a new enemy submarine is in existence. The simulation can be used to 
indicate how well the present vehicles will perform against the new 


threat relative to the old threat. 
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Je Conclusions. 

The example in Section 4 is the type of comparison that can be per- 
formed for various kinds of ASW vehicles. There are several decisions-- 
evasion tactics used, submarine type used, and so on--that have to be 
made prior to using the program; but once the decisions are made for a 
given situation, they have to remain the same for all systems being com- 
pared. Several modifications and suggestions for the model are discussed 
below. 

It might be argued that an elliptical coverage function should be 
used instead of a circular coverage function because of the elliptical 
sensor locating error. This change can be made by rewriting SUBROUTINE 
OCIP of the computer program (Appendix C). 

Instead of evasion turns Limited to 90 degrees either side of the 
course line, they could be programmed for 360 degrees evasion; i.e., 180 
degrees either side of the course line. 

It must be realized that the depth of the submarine has no direct 
influence on this simulation. Therefore, this has to be kept in mind 
when parameter values such as EXPRNG are put into the simulation. It is 
possible that EXPRNG could be less in an evasion problem where the sub- 
marine is more likely to increase its depth, and therefore, EXPRNG should 
be reduced. 

It would be possible to give the vehicle freedom of movement before 
transit commences, rather than being motionless. But the advantages of 
the passive tactic for the vehicle are then relinquished. 

The restriction on Tl in the alerted submarine case could be abol- 
ished by programming for the possibility of having a speed change exe- 


cuted before a turn is completed. 
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Possibly a passive sensor detached from the vehicle could transmit 
current information about the submarine's actions while the vehicle is 
inbound to intercept the submarine. This would increase the probability 
of success. 

In the foregoing sections, the simulation has been described and the 
mathematical development shown. In addition, the use of the computer 
program has been illustrated and an example and its results briefly ana- 
lyzed. The use of the simulation is simple and the results easily under- 
stood. Therefore, it is recommended that the use of this simulation be 
employed to compare and analyze proposed ASW vehicles and operational 


ASW vehicles. 
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Notation 


A(1),.-.,A(Q5) 


AB 


BETA 


BIAS(N) 


BRGER 


CLDT 


DELTA 


DIST 


EP 


EXPRNG 


NEXT 


NU 


PAVG 


PRCNTN 


APPENDIX A . 
GLOSSARY 
Meaning 


Number of standard deviations from the mean of 
the normal density function. Each A(i) is the 
midpoint of a segment of area which represents 
20% of the total area under the normal density 
function. The A(i)'s are used in the predic- 
tion error of the unalerted submarine case. 


Evasion turn in degrees from the estimated 
submarine course. 


True bearing of the submarine from the vehicle 
(degrees). 


Distance from EP to the terminal point of the 
nth evasion. 


Sensor bearing error in degrees. 

Land time (min.) which includes time to land 
into the wind and taxi time to get into posi- 
tion to employ active sensors. 

Bound to the right of the known submarine 
course. This bounds to the right the set of 


possible evasion turns (degrees). 


Distance the submarine travels during blind 
time, assuming a constant course and speed. 


Estimate position of the submarine after blind 
time has expired. 


Expected range of the active sensor used by the 
vehicle in re-establishing contact at the EP 
(in nautical miles). 


Number of data runs for the computer program. 


Number of increments considered with the 
(PSI-DELTA).range of the set of possible turns. 


Average probability of success for N evasions 
considered in the alerted submarine case. 


Navigation error in percent of distance trav- 
eled. 
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Notation 
PRCNTS 


PROB 


PSI 


PTOT 


RELT 


RT 


SBIAS(I,J) 


SIGJT 


TB 


TD 


Tl 


THDEL 


THETA 


TOT 


Meaning 
sensor range error in percent of sensor range, 


Probability of success for the unalerted sub- 
marine case. 


Bound to the left of the last known submarine 
course. This bounds to the left the set of pos- 
sible evasion turns (degrees). 


Pre take-off time (min.) which commences when 
the vehicle can no longer observe the submarine 
and includes decision time, warm-up time, and 
sensor retrieval time. 


Range from the vehicle to the submarine at the 
last known position (nautical miles). 


Target relocation time (min) which includes the 
time to deploy sensors and ends when the sub- 


marine is relocated. 


Reaction time which is the sum of PTOT, TOT, 
CLDT, and RELT. 


Distance between the ci, 4) th cell and EP. 
Total standard deviation of errors. 

Blind time which is the sum of the reaction 
time (RT) and the transit time (TD). 


Time it takes the vehicle to travel from its 
initial position to EP. 


Tactical interval of time (min.) after evasion 
commences during which the submarine travels at 
fast speed. After Tl has elapsed, the sub- 
marine slows to silent speed. 


Submarine course error (degrees). 


Course of the submarine (degrees) at the last 
known position. 


Take-off time (min.) which includes take-off 
roll.and the time to maneuver into the wind. 


Turning radius of the submarine. 
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Notation 
UERR 


UNU(1),..-,UNU(NU+ 1) 


Ul 
U2 
U3 
V 


WIPAVG 


(X2,Y2) 


(X(N), Y(N)) 


Meaning 

Submarine speed error (knots). 

The set of probabilities for N=l1,..., NU+l 
evasion turns. These are used when evasion 
turns considered are not equally likely. 
Submarine cruise speed (knots). 

Submarine silent speed (knots). 

Submarine fast speed (knots). 


Average vehicle transit speed (knots) 


Weighted probability of success for the alerted 
submarine case. 


Last known position of the submarine. 


Terminal point of the nth evasion after blind 
time has expired. 


Average vehicle transit speed (knots). Neglect 
acceleration and deceleration time. 
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APPENDIX B 


INPUT AND 
UNIT CHANGES 





INPUT 
PARAMETERS 









CHANGE SPEEDS 
FROM KNOTS 
TO 


NAUTICAL MILES 
PER MINUTE 


CHANGE ANGLES 
FROM DEGREES 
TO 
RADIANS 











MANEUVERING BOARD 


© PROBLEM 


Xl= Ul x SIN (THETA) 
Viz Ul ¥ COS(THETA) 


X2=RAI x SIN(BETA) 
Y2= RAT x COS(BETA) 


E a WYXEY V2 - (x0)? 
X27 0 . 
E (Ro) Y32)= -Y¥3(1) 
Or) ABIYY = ]¥2 -Y3(1)] 


AB2YY= |Y2-Y3(2)] 






Y2=O 


T 


X3(1) =] v2 -(¥i)* 


X3(2)= — X3(1) 


: 70 =O 
ABIXX= }x2 - x30) ABIYY— AB2YY: 


AGIXX=]¥2 —X3(2)| <0 








Le=2 


5 XA(L) =X} 
=O f anix-Aap2xx )22 


<0 


ERROR 
STATEMENT 


ERROR 
STATEMENT 


NO=NO+ | 
NM=Ms] 





NO=NO + | 
M=M+| 





Is] 
y30)=¥1 (§) 
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_-2 xXx (Y2)4 eV (¥ 
Aa + (¥2)* Ax2)* 
ae we _2xf JENIN + (vi)? -y? 












_-B+Bt-4AC- 


X3{)= -yN 


~B-Y¥B*-4AC 
30 2A 






ABSXK |= |¥2-X3(1)] 
| ABsKx2 = |X2-x3(2)I 


*40 >oa 


Note oval eu Y3¢= -Y33 


ERROR 
STATEMENT 








agsyv |= |Y2 -¥33] 


AGSYY \={| Y2~ Y3) 
| | ABSYY2 = | ¥2-Y34| 


AGSYY2=| Y2-Y32] 





NO=NO+\ 
M=M+\ 


re. ae ABSYV) — ABSYY2 
Aree ne 
ERROR 


ERROR 


STATEMENT STATEMENT 









Lz2 





I= | 
Y3(L)=Y32 


¥3(I) =Y34 









Te2 
¥3(1) = Y33 


o. 
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r= 
¥3(I)=Y3\ 











RELSPOD = 300) -X 1)" + (¥3(2)- Y1)> 





RELATIVE 
SPEED 


TD = RAT /RELSPD 


RT = PTOT + TOT 
+CLOT + RELT 


IX=0 
SUBT= RAT /U3 


REACTION TIME TOD 
GREAT. SVG CAN 
KILL VEHICLE. 






TR=TO+r aT 
DIST = V3xXTB 


X4 = DIST x SIN (THETA) 
+X2 


Y4 = DIST x COS (THETA) 
+Y¥2 


DISDAT=/(xy)+(¥4)* 


BLIND TIME 
OISDAT 
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ALERTED CASE 







ae oh ERROR 
as STATEMENT 
q0°~ PSI ° 
x x 
ae NO=NO+] 
M=M+l 


ETA= 


SIGNAV=PRONTW X RAT 
STGSR=PRCNTS x RAT 
GSBR=RAI x SIN(BRGER) 


SIGAPR= VSIGSA xSIGSBR 





>O ; 
TE- 
£0 
<0 
ory ERROR 
STATEMENT 
>0 
=| 


: IFLAG 


Th=TI 


Paes 











UN= NU +1 
NUN= UN 





DN= | DELTAL -3¢0"( 


£0 
DN = .0009 
ELTAL = 360° : 
al DELTA(N)=0 
TTT =O 
A&=@Qe0 
STR=0 





.e) 







70 


6 
i DELTAL - 260 
DELTAL = DELTAL-360° 


I 
i 


i Eo) 


OM= [DELTAL -90°| 





Z 
4,79 
<Q = 6 +0 
SEEDY pean 
ii -|=0 DELTA(N)= AB 
DELTA(N) = DELTAL . ar, STRXAB 
AB= BB = DELTAL AG=DELTA(N) = 90° = 3 
X(4} = STR 
YN) =STR+(TB- 
EVASION TIT)* U3 EVASION 





TURN | TURN 


EVASTON 
TURN 
(2) ¢ v 











X(N) =STRX SIN (8B) +U3x(TB-TTT)x Cos (AB) 
Y(N) = STR=STR x COS (QB) + U3 x (TE-TTT) x SIN (AB) 





XXX = X(N 
ERROR nt IFLAG 22 YYY = in 
STATEMENT | 





TR= SAVE -TI 


WO=NO + | 


MeM+| X(N) © TB x U2 x CoS (AB) 


YN) =18.x U2 x SIN (AB) 


(=) - X(N) =XXX + X(N) 


Y(N)2NVY + YOM) 


COMPUTE PPK(N) 
FOR ALERTED CASE 












XX = 
YY = 


OTST 
6) 


BIAS (N) =) K(N)- orst)\+ (¥(W))™ 





O= nras(N) J Sica 
AKTLL = EXPRNG /SIGIT 


COMPUTE EVASION PROGAGILITY OF SUCCESS , Pk, 
FOR THE NG EVASION VSING THE CIRCULAR 
COVERAGE FUNCTION, 





PPK(N) = Pk 
Sum = sum + PPIK(N) 
SUMM = SUMM + PPK (N) x UNY (N) 


PROBABILIES 
OF S¥CccESS 
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WTPAVG ® sun /ux 


PAVG = sum / UN 


ry 


AVERAGE EVASION 


SUCCESS PROB. 





OVUMUNU = O. 


SUMUNU =SUMUNU +UNU(N) 


| 
4 
4 


WEIGHTED EVASTON 70 
success eros. 


SUMUNU 


TT 
© 


ERROR 
STATEMENT 





NO=NO+ | 
A= M+] 
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“WEIGHTED PROG, 
Nor USED ™ 











UNALERTED CASE 


W= (UI -AA()xUERR) x TB 


THETER= AA(S) x THDEL 


% 


> 


SUBX(1,3)=VV x COS (THETER) 
SuBY (1,5)=VV x SIN(THETER) 


BIAS(153)=\/60(—- SuGX(I,)\74 (VY- SUBY(1,3))” 


D= BIAS(I,3)/SIGIT 
OAKILL = EXPRNG/SIGIT 


. 
“ 


COMPUTE UNALERTED PROBABILITY OF SUCCESS, PK FOR (1,J)T# 
CELL USING CIRCVLAR COVERAGE FUNCTION. 


SUM <= SUM + PK 


PROG = SUM X .04 


UNALERT ED 


i (400) 


SucCCESS PROG. 
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SUBROUTINE INPUT 


SUM = O 
SUMM= 0 


READ 26 INPUT 
PARAMETERS 





ISK 
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